Consistent cosmological modifications to the Einstein equations 
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General Relativity (GR) is a phenomenologically successful theory that rests on firm foundations, 
but has not been tested on cosmological scales. The advent of dark energy (and possibly even the 
requirement of cold dark matter), has increased the need for testing modifications to GR, as the 
inference of such otherwise undetected fluids, depends crucially on the theory of gravity. In this work 
I outline a general scheme for constructing consistent and covariant modifications to the Einstein 
equations. This framework is such that there is a clear connection between the modification and 
the underlying field content that produces it. I conclude by a simple metric based modification 
of the fluctuation equations for which the background is exact ACDM and present its impact on 
observations of the cosmic microwave background radiation. 



The theory of gravity plays a fundamental role in our 
modelling and understanding of the universe. If we are 
to know the matter constituents of the universe, we 
have to be sure we understand what is the underlying 
gravitational theory. Einstein's General Relativity (GR) 
has played a key role in formulating modern cosmology, 
first as a smooth Friedmann-Lemaitre-Roberson- Walker 
(FLRW) spacetime, then at the level of linearized fluctu- 
ations about this spacetime. 

General relativity is a very solid "principle" theory 
from the theoretical point of view, (and quite understand- 
ably the aesthetical point of view). The Lovelock-Grigore 
theorem |l|, |2| asserts that GR with a cosmological con- 
stant is unique under the following assumptions : geom- 
etry is Riemannian and the gravitational action depends 
only on the metric, it is local and diffeomorphism invari- 
ant and leads to 2nd-order field equations. Relaxing any 
of these assumptions can lead to more general gravita- 
tional theories, e.g. adding extra fields [3, 0,0, having 
higher derivatives [7|, having a pregeometry [8§ or making 
the theory non- local [9J- This is not an exhaustive list but 
possible theories fall into one or more categories above. 

However as nice as we may think that GR is, the ul- 
timate judge is experiment. Indeed, different aspects of 
GR have been vigorously tested in the lab, in the solar 
system and with binary pulsars, all of which lie in the 
strong curvature regime (compared to cosmology). 

The discovery that the expansion of the universe is 
accelerating opens the possibility that general relativity 
breaks down on large scales or low curvatures. It may also 
be that the apparent missing mass in the universe is not 
in the form of a cold dark matter particle but once again 
due to departures from general relativity. This opens the 
need for cosmological tests of gravity, and much work has 
been carried out in this direction [T(| HlxI. Il2| at various 
levels. More recently Hu and Sawicki |13j . and Hu [Til ], 
have laid down a fully covariant formulation of modifica- 
tions to gravity under well motivated assumtions. 

In this work I outline a general scheme for constructing 



consistent modifications to the Einstein equations. The 
scheme is such that, one can clearly classify the modi- 
fications according to whether they obey or violate dif- 
feomorphism invariance, need extra fields, or stem from 
higher derivative theories. Indeed the advantage of this 
method is the direct connection between the field con- 
tent and the modifications. Additional assumptions as 
in [3, [13, [ll| can always be used at the very end but we 
shall not consider this possibility here. I close the paper 
by constructing the most general modification for which 
the FLRW background is exact ACDM and illustrate the 
effects on obscrvables in a simple subcase. 

Following Hu [14[ we start by putting the gravitational 
field equations in the form 



Gab = 8ttGT 



U ab 



(1) 



Here G ao is the Einstein tensor for the universal mat- 
ter metric g a b, X^™ " 1 "' is the stress-energy tensor of all 
known forms of matter (like baryons, photons and neu- 
trinos), U a b is a general tensor which encapsulates all the 
unknown fields/modifications and can depend on T a b for 
each field and various combinations of metric functions 
(such as curvature tensors). The assumption of upto 2nd 
order field equations, translates to having only upto first 
derivatives of the extra field in Uq v and upto 2nd deriva- 
tives in Uij . Relaxing this assumption is possible and will 
simply give higher order field equations, but one has to 
be cautious that quite generally higher derivative theories 
lead to instabilites. 

At this point we have to decide on the field content, i.e. 
whether U a b depends on additional fields or the metric 
alone. In the former case, we must add one field at a 
time and ensure its energy conservation by applying the 
Bianchi identity which directly translates to 



V a U a b = 



(2) 
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and gives a field equation for the extra field. Violating the 
Bianchi identity leads to entirely arbitrary parametriza- 
tions and will not be considered. The caveat of this ap- 
proach is that in the case that more than two indepen- 
dent degrees of freedom are present in U a b, one would 
have to supply extra equations for these, not given by 



2 



the Bianchi identity (see [H, Oj| [l6[ for examples) . Fi- 
nally, there could be interactions between fields in T a ^ 
and U a b- For simplicity I do not consider this possibility 
further, but it is straightforward to add. 

We now split the dynamics of the problem in the back- 
ground FLRW dynamics and their fluctuations about 
that background. The FLRW metric is ds 2 = a 2 (—dr 2 + 
qijdx' l dx J ) where r is the conformal time, a is the scale 
factor and qij is a spacial metric of constant (dimension- 
less) curvature K. The FLRW assumption means that 
effectively we are considering a collection of scalar fields 
on the spacial hypersurface of homogeneity and isotropy. 
This boils down to requiring that the Lie derivative of 
the extra field vanishes for all six Killing vectors of the 
FLRW spacetime. Examples are a scalar field 4>(t), a 
vector field with components = (A(t), 0,0,0) and a 
tensor field X a b whose only non-vanishing components 
must be X° = -X(t) and X l j = Y(t)S i y In the special 
case of a unit-timelike vector field A(t) is pure gauge and 
the contribution of such a field to the FLRW equations 
is generally given in terms of functions of a, a and d Q. 

Lets define Ep = — a 2 G° and Er such that a 2 G l ■ = 

■ 2 

Er8 % j, which are explicitely given by Ep = 3^ + 3K and 

. 2 

Er = + The FLRW equations corresponding 

to (JTJ> are then simply written as Ep = &nGa 2 J^i Pi + 
a 2 X and E R = %-kGo 2 J2 t Pi + a 2 Y. Applying the 
Bianchi identity then gives Ep + ^{Ep + 3Er) = 0. and 
the fluid equation p + 3-(p + P) = 0. Applying the 
Bianchi identity on U a b imposes 



X - 
V 



8 - 

E e 
E P 



•X+k [6f£ + 2V 2 V>] 
^E P + l\E R -2iE R \Z 



w -^w + 



a a 
1 



6^ + 6(|-2^ 



n 



+2V 2 £ 

♦ n + i' 



w - 3w(l+w)i] £ 



E a ^Ea-~Z(E f + E h )£ 



TABLE I: Gauge transformations for the metric, fluid and 
Einstein tensor variables. 



perturbation II such that 6T*j = HpS 1 j and shear S, 
such that the shear tensor is = DyE. 

We now decide whether the parametrization should 
obey diffeomorphism invariance. At the linearized level 
this is the requirement that all field equations must be 
gauge form-invariant. Let us demonstrate what is gauge 
form-invariance for the case of standard GR coupled to 
a fluid. Gauge transformations are infinitesimal diffco- 
morphisms generated by a vector field £ a which can be 
parametrized as £ M = o(— £, Viip). All perturbations 
apart from E above are not gauge invariant but transform 
with £ and ip as in table I. Consider the 8G° Einstein 
equation which is 



i(V 2 + 3iq 



X-V 2 v 



X + 2V 2 C 



X + 3-{X + Y) 
a 







(3) 



The above equation will give the background equation for 
the extra field, or additional constraints on X and Y in 
the absence thereof. We see that any FLRW background 
can be modelled via one arbitrary function Y(t) and a 
second function X(t) found by solving ([3]). 

Adding a scalar field amounts to letting X = 
X((/>,<f>,a,a) and Y = Y(<f>, <fi, <f>, a, a, a), if we are to ex- 
pect at most 2nd order field equations for the scalar. 
For example a canonical scalar field corresponds to X — 
2±2-0 2 + V(cj)) and Y = ^<j) 2 - V(<p) while the Jordan- 
Fierz-Brans-Dicke theory [3] to X = 87rGa 2 (e* -l)p + 
" 12 + M4> and Y = ^Ga 2 (e* - 1)P _ + i0 - 
2 ). 

Consider now the (scalar) fluctuations about the 
FLRW metric as 



2a 2 " 
2 



ds 2 = -a 2 (l - 2E)dt 2 - 2a 2 (Vi()dtdx l 
+a 2 [{l + -x)q%j + D ij u]dx i dx j 



(4) 



where D.y = V^Vj — |ffy V 2 is a spacial traceless deriva- 
tive operator. A perfect fluid is described at the fluctua- 
tion level by a density contrast 8 , momentum 6 such that 
its total momentum is it, = aV \6, dimensionless pressure 



-6— E = 8irGa 2 p6 



(5) 



If we perform a gauge transformation X — » X' to all 
perturbations above, the SG° Einstein equation becomes 



l(V 2 + 3K) [x'-VVl +- 
3 L Jo 

+6^S' = 87rGa 2 p5' 



id + 2V 2 C' 



(6) 



i.e. it retains its exact form (it is form-invariant), with 
the only change being a relabelling of the perturbations 
with X — > X' for each variable. Gauge form-invariance 
always holds for all field equations which stem from a 
diffeomorphism invariant action, no matter how compli- 
cated the theory is. 

We can shortcut testing for gauge form-invariance as 
follows. First define the three gauge non-invariant po- 
tentials y = z> + 20 J = x- V 2 v and W = \ + 2V 2 C = 
J + W 2 V, which are the only three combinations of met- 
ric variables appearing in the perturbed Einstein tensor. 
They transform only with the gauge variable £. Then de- 
fine the two gauge- invariant potentials <3E> = — \J + \ 

and i> = —E - \V - . We can now split the per- 

turbed Einstein tensor into a gauge invariant and a gauge 
non-invariant part which involves the variable V . For 
simplicity let us define Ea = —a 2 8G° Q , Eq is such that 
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-a 2 SG° i = ViE @ , E P = a 2 6G\ and E s is such that 
a 2 [5G l } - \§G k k 8 i j \ = D^Es. Explicitely we get 

E A = 2(V 2 + 3A)$- 6- ($ + -*) 

a a 

A±(E F + E R )V, (7) 
2 a 

£ e = 2(6+-*) + J(B F +^Ji)V; (8) 
a 2 

Ep = e4- + +12-4 + -*) 

at \ a J a a 
-2(V 2 + 3iT)(<£> - *) - 3(£ F + E R )i< 

+ l[E R -2^E R ]V, (9) 

and 

£ s = l> - *. (10) 

The Einstein equations are then given as 

E A = 8TrGa 2 pS + U A (11) 

E e = %-KGa 2 {p + P)6 + U@ (12) 

E P = 2AnGa 2 pn + U P (13) 

£ s = 87rG(p + P)S + C7 E (14) 

where the Ui variables are defined in the same way as for 
the Ei variables with SG a b replaced by 5U a b . 

The importance of gauge form-invariance cannot be 
overemphasized. In fact gauge form-invariance severely 
constrains the terms involved. For had we decided to let 
the RHS of the Einstein equations to be SirGa 2 fi(r)S, 
87rGa 2 / 2 (r)0 and 247rGa 2 /3(r)II we would have found 
that by virtue of the background equations, the Einstein 
equations are gauge form-invariant if and only if j\ = 
/.3 = p, and fi = p + P. A less trivial example is when 
X = —\ for some scale r c . If we decide that this term is 

due to metric modifications, then U A = — (i J+-3) and 

^ e = — str" (jj "t" a")' If on the other hand we posit that 
it is due to a unit-timelike vector field perturbed as ^4^ = 
aSliOi then a possibility is U A — — 2^)a. In the 
last case we pick higher derivatives, and a gauge- invariant 
counter-term must be added to cancel them |17| . 

Once we have added the gauge non-invariant terms and 
correctly fixed the functions multiplying them by requir- 
ing gauge form-invariance to hold, we can proceed to add 
more gauge-invariant terms involving the extra fields and 
the gauge invariant potentials $ and *. If we want to 
consider only parametrizations which lead to 2nd order 
field equations in all variables, then we have to be care- 
ful what terms we add and where. For example, when 
we add metric terms, we can add upto first derivatives 
in the two Einstein constraint equations (jTT|) and (fl"2")) . 
and up to second derivatives in the two Einstein propa- 
gation equations (fTB")) and (fT4")) . Since $ is of first order 



in the metric variables while * is of second order, we can 

add $ in all four Einstein equations but $ and * only in 
the two propagation equations (TIB")) and (fl"4"|) . Note that 

although $ + 2 * is of lst-order in the metric perturba- 
tions, it contains 2nd derivatives of the scale factor, and 
so it cannot be added to (fTT|) and (Tl2|) . unless we relax 
the 2nd-order field equations constraint. The reason it 
is allowed in the definition of E A above is because, there 
are also 2nd derivatives of the scale factor appearing in 
the gauge non-invariant part of E A proportional to V. 
Thus when E A is written in terms of the actual metric 
potentials, the a terms cancel and the final expression 
contains only first derivatives in all the variables. On the 
contrary, when we add solely gauge-invariant terms we 
no longer have this luxury. 

Let us also note that the gauge variable ip is not in- 
volved in the transformation of the Einstein tensor. Thus 
if any extra field does transform with ip it will always ap- 
pear in combination with u, £ or \ m the field equations, 
in a way that the whole combination does not transform 
with tp. An explicit example can be found in (l6j . 

We finally utilize the Bianchi identity which at the 
linearized level simply translates in terms of the added 
variables U A , U&, Up and Us as 

U A + -U A -V 2 Ue + la 2 (X + Y)W+-U P = (15) 
a 2 a 

and 

U e + 2-U e - \u P - \{V 2 + 3K)Us 
a 6 3 

+a 2 (X + Y)E = (16) 

This provides us with the field equations for the extra 
fields [l3|, or with additional constraints on the added 
functions in the absence thereof. 

I now illustrate the above scheme, by finding the most 
general diffeomorphism invariant modification to the Ein- 
stein equations for which the background cosmology is 
the plain ACDM model, no extra fields are present, 
and no higher derivative than two is present in the field 
equations. Since there are no extra fields and the back- 
ground is unchanged from ACDM we can only add 
gauge invariant terms to Einstein equations by setting 

U A = A$, U e = Up = Ci* + C 2 <l + C 3 * and 

Ut, = £>i* + £> 2 $ + Z>3*, for operators A, B, d and 
T>i. Applying the Bianchi identity we get two equations 

involving $ and * and consistency requires that these 

equations must be satisfied whatever the values of <t, <& 
and *. A sufficient condition is found by setting the co- 
efficients of these terms to zero which gives C3 = 2?3 = 0, 
A = -fC 2 , B = \C 2 + f (V 2 + 3A')£> 2 and the two 
differential equations A + %A - V 2 B + |Ci = and 
B + 2fB - ±Ci - |(V 2 + 3A)2?! = 0. A quick exam- 
ination reveals that if A and B are both zero then we 
get exact GR. The same holds for T>\ and P 2 , hence a 
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FIG. 1: Upper panel: The CMB spectrum for the simple mod- 
ified gravity model in the text. The solid curve is the plain 
ACMD model ((3 = 0), while the dotted, dashed and dot- 
dashed curves are with j3 = {0.1, 0.5, 1} respectively. Lower 
panel : The time evolution of g — (as in [13, (HI) at 

k = 10~ 3 A/pc _1 for the same set of models. 

generic prediction of this kind of modification to GR is 
that $ — ^ should deviate from the GR value. 

Lets illustrate the effect on observables in a simple sub- 
case for which B = C x = V x = 0, A = SSL, C 2 = 

3H 2 l 

and 2?2 = ~2^"^- Thus, we parametrize deviations from 
GR with a single dimensionless parameter (3. We let 
7 = 2k 2 a+i3H' i anc ^ ml< ^ ^ e Einstein equations in the 
synchronous gauge (defined as S = Q = and X = h, 



—k 2 v = h + 677) . The h equation is 

-h= (1 - 1 )8nGa 2 pS + 2k 2 r 1 -6 1 -r] (17) 
a a 

while the t) equation is unchanged. 

The perturbation equations were solved numerically in 
both the synchronous and in the conformal Newtonian 
gauge for numerical consistency. The upper panel of Fig. 
1 shows the CMB angular power spectrum 1(1 + l)Ci for 
a ACDM universe (j3 — 0) contrasted with non-zero (3. 
We see that for this particular model, the effect of non- 
zero P is to decrease power on large scales including even 
the first peak. The lower panel of Fig. 1 shows the 
time variation of q = f— i for the same set of models at 

k = 10 _3 Afpc _1 . We see that like other modifications to 
gravity, the effect is to make g grow. In contrast to other 
parametrizations of modified gravity [l2l.[l8| however, the 

difference of $ — ^ is sourced by <SE> rather than 

This scheme can be used to consistently add terms to 
the Einstein equations which may or may not depend on 
extra fields, and possible metric functions. The fluctu- 
ation equations can then be solved to obtain observable 
spectra. An important issue is to find predictions con- 
cerning these extra parameters for upcoming precision 
experiments such as the Square Kilometer Array (l9| and 
the EUCLID project. 
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